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ABSTRACT 

We report the results of interferometric HCN(l-O) and HCO + (1-0) observations of four luminous 
infrared galaxies (LIRGs), NGC 2623, Mrk 266, Arp 193, and NGC 1377, as a final sample of our 
systematic survey using the Nobeyama Millimeter Array. Our survey contains the most systematic 
interferometric, spatially-resolved, simultaneous HCN(l-O) and HCO + (1-0) observations of LIRGs. 
Ground-based infrared spectra of these LIRGs are also presented to elucidate the nature of the energy 
sources at the nuclei. We derive the HCN(1-0)/HCO + (1-0) brightness-temperature ratios of these 
LIRGs and confirm the previously discovered trend that LIRG nuclei with luminous buried AGN 
signatures in infrared spectra tend to show high HCN(l-O) /HCO + (1-0) brightness-temperature ratios, 
as seen in AGNs, while starburst-classified LIRG nuclei in infrared spectra display small ratios, as 
observed in starburst-dominated galaxies. Our new results further support the argument that the 
HCN(1-0)/HCO + (1-0) brightness-temperature ratio can be used to observationally separate AGN- 
important and starburst-dominant galaxy nuclei. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: ISM — radio lines: galaxies - 
galaxies: individual (NGC 2623, Mrk 266, Arp 193, and NGC 1377) 



1. INTRODUCTION 

Luminous infrared galaxies (LIRGs) radiate the bulk 
of their large luminosities (L > lO 11 ^©) as infrared dust 
emission. Large infrared luminosities mean that: (1) lu- 
minous energy sources are present, but are hidden be- 
hind dust; (2) energetic radiation from the hidden energy 
sources is absorbed by the surrounding dust; and (3) the 
heated dust grains re-emit this energy as infrared ther- 
mal radiation. To understand the nature of LIRGs, it is 
essential to unveil the hidden energy sources; namely, to 
distinguish whether starbursts (= nuclear fusion inside 
stars) are dominant, or active galactic nuclei (AGNs; ac- 
tive mass accretion onto a central compact supermassive 
black hole [SMBH] with >10 6 M Q ) are also energetically 
important. 

Unlike AGNs surrounded by toroidally shaped dust, 
which are classified optically as Seyferts (Veilleux & Os- 
terbrock 1987), concentrated molecular gas and dust in 
LIRG nuclei (Sanders & Mirabel 1996) can easily bury (= 
obscure in virtually all directions) the putative compact 
AGNs, making AGN signatures very difficult to detect 
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optically. However, it is crucial to understand the ener- 
getic role of such optically elusive buried AGNs in LIRG 
nuclei. 

A starburst (nuclear fusion) and buried AGN (mass 
accretion onto a SMBH) have very different energy gen- 
eration mechanisms. Specifically, while UV emission is 
predominant in a starburst, an AGN emits strong X- 
ray emission in addition to UV. Additionally, in a nor- 
mal starburst, the stellar energy sources and dust are 
spatially well mixed, so energy sources are spatially ex- 
tended. Since the energy generation efficiency of the nu- 
clear fusion reaction is only ~0.5% of Mc 2 (M is the mass 
of material used in the nuclear fusion reaction), the emis- 
sion surface brightness of a starburst is modest and has 
both observational (Soifer et al. 2000) and theoretical 
(Thompson ct al. 2005) upper limits (~ 10 13 L Q kpc~ 2 ). 
On the other hand, in an AGN, the mass accreting SMBH 
is spatially very compact, and thus more centrally con- 
centrated than the surrounding gas and dust. The high 
energy generation efficiency of an AGN (6-42% of Mc 2 ; 
M is the mass of accreting material; Thorne 1974) can 
generate large luminosities from a very compact region, 
producing a very high emission surface brightness (Soifer 
et al. 2000). These differences between an AGN and a 
starburst could create differences in the properties of the 
surrounding molecular gas and dust that may be distin- 
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guishable based on observations at wavelengths of low 
dust extinction. 

Molecular gas emission lines in the millimeter wave- 
length range can be used effectively to investigate the 
hidden energy sources of LIRG nuclei. First, dust ex- 
tinction is very small. Second, theoretical calculations 
predict that X-ray- and UV-emitting energy sources show 
different behaviors of molecular line emission due to dif- 
ferent chemical reactions (Mcijcrink ct al. 2006). Fi- 
nally, mid-infrared 10-20 /xm emission is strong in an 
AGN because of hot dust thermal emission heated by 
the high surface brightness energy source of an AGN. 
The strong mid-infrared emission could selectively en- 
hance particular molecular line emission through an in- 
frared radiative pumping mechanism (Aalto et al. 1995; 
Garcia-Burillo et al. 2006; Guelin et al. 2007; Weiss et 
al. 2007; Aalto et al. 2007). In this paper, we focus 
on HCN(l-O) and HCO+(1-0) lines, because an AGN 
could enhance HCN(l-O) emission through an increased 
HCN abundance (Lintott & Viti 2006) and/or infrared 
radiative pumping (Aalto et al. 1995; Garcia-Burillo ct 
al. 2006; Weiss et al. 2007). In fact, Kohno (2005) found 
observationally that HCN (1-0) emission is stronger, rel- 
ative to HCO + (1-0), in AGN-dominated nuclei than in 
starburst galaxies, demonstrating the potential of this 
HCN(l-O) and HCO+(1-0) based method for the pur- 
pose of distinguishing between an AGN and a starburst. 
This method may be effective even for a Compton thick 
(N H > 10 24 cm~ 2 ) buried AGN, the detection of which 
is very difficult with direct X-ray observations. 

HCN(l-O) and HCO+(1-0) observations of a large 
number of LIRGs have been made, using single-dish ra- 
dio telescopes (Gao & Solomon 2004; Gracia-Carpio et al. 
2008; Krips et al. 2008). However, the beam sizes of these 
observations are so large (>25 arcsec) that severe con- 
tamination from spatially-extended star-forming emis- 
sion is unavoidable, possibly hampering the detection 
of the AGN signatures of LIRG nuclei, where putative 
AGNs are expected to be present. Furthermore, LIRGs 
often show multiple-nucleus morphologies with small sep- 
arations (<10 arcsec), for which spatially-resolved study 
is not possible with single-dish radio telescopes. Fi- 
nally, the HCN(l-O) and HCO+(1-0) data were col- 
lected at different times under different weather condi- 
tions, possibly increasing uncertainty regarding their rel- 
ative strengths due to the ambiguity of inter-calibration 
of data. 

To overcome these issues, we performed millimeter in- 
ter ferometric, spatially-resolved, simultaneous HCN(1- 
0) and HCO + (1-0) observations of LIRGs, using the 
Nobeyama Millimeter Array (NMA) (Imanishi et al. 
2004, 2006b; Imanishi & Nakanishi 2006; Imanishi ct 
al. 2007b). In our study, only nuclear HCN(l-O) and 
HCO + (1-0) emission was extracted to investigate their 
ratios in brightness temperature (oc flux x A 2 ). The 
sample contained mostly ultraluminous infrared galax- 
ies (ULIRGs; L m > 10 12 L Q ; Sanders et al. 1988a) 
showing luminous buried AGN signatures at other wave- 
lengths, and we found that their HCN(1-0)/HCO+(1 
0) brightness-temperature ratios were generally higher 
than starburst galaxies with Lir < 10 12 L©. This re- 
sult supports the presence of luminous buried AGNs 
in these observed ULIRGs. However, scenarios that 
try to explain large HCN(1-0)/HCO+(1-0) brightness- 



temperatures ratios, without invoking luminous AGNs, 
have been present. Examples include the high turbu- 
lence (Papadopoulos 2007) or high density (Imanishi 
et al. 2007b) of molecular gas. In fact, ULIRGs have 
very dense, highly disturbed molecular gas compared to 
normal LIRGs with Lir < 10 12 L Q (i.e., non- ULIRGs) 
(Sanders & Mirabel 1996; Soifer et al. 2000, 2001). 
Observations of normal LIRGs (non- ULIRGs) , with or 
without buried AGN signatures, would help elucidate 
whether the high HCN(1-0)/HCO+(1-0) brightness- 
temperature ratios in observed ULIRGs are indeed 
caused by luminous buried AGNs, or are simply due to 
different molecular gas properties in ULIRGs. 

In this paper, we present the results of interferometric, 
simultaneous HCN(l-O) and HCO + (1-0) observations of 
such normal LIRGs (non-ULIRGs), to investigate the 
HCN(1-0)/HCO+(1-0) brightness-temperature ratios in 
the nuclear regions. This is the final paper of our series 
of interferometric HCN(l-O) and HCO+(1-0) surveys of 
LIRGs using NMA (Imanishi et al. 2004, 2006b; Iman- 
ishi & Nakanishi 2006; Imanishi ct al. 2007b). Ancillary 
ground-based infrared K- (1.9-2.5 /im) and L- (2.8-4.1 
fim) band spectra are also utilized to better understand 
the hidden energy sources of these millimeter-observed 
LIRG nuclei. Throughout this paper, we adopt H = 75 
km s _1 Mpc -1 , Qm = 0.3, and f^A = 0.7, to be consistent 
with our previously published papers. 

2. TARGETS 

We observed the four LIRGs, NGC 2623, Mrk 266, Arp 
193, and NGC 1377, primarily because they were non- 
ULIRGs with or without buried AGN signatures and so 
were appropriate sources to address the above-mentioned 
issue (§1), and because we estimated that we could de- 
tect HCN(l-O) and HCO+(1-0) emission with NMA in 
a reasonable amount of telescope time. Table 1 sum- 
marizes the infrared emission properties of these LIRGs. 
An angular scale of 1" corresponds to a physical size of 
0.12-0.53 kpc at the redshifts of these LIRGs (z = 0.006- 
0.028). 

NGC 2623 (z = 0.018) is a LIRG (L m = 1O 115 L ) with 
two prominent merging tails along the north-eastern and 
south-western directions (Rothberg & Joseph 2004). It 
is optically unclassified (Keel 1984; Veilleux et al. 1995); 
no obvious AGN signature exists in the optical spectrum. 
Its infrared 5-35 /zm spectrum obtained with Spitzer IRS 
is dominated by polycyclic aromatic hydrocarbon (PAH) 
emission, as is usually observed in starburst galaxies 
(Brandl et al. 2006). However, high spatial resolution 
ground-based infrared 10-/xm imaging observations re- 
veal that the bulk of the 10-/im emission in NGC 2623 
comes from compact nuclear cores with <400 pc, the 
surface brightness of which can be as high as ^10 13 L Q 
(Soifer et al. 2001). This value is close to the maximum 
found in starburst phenomena (Soifer et al. 2000). The 
2-10 keV X-ray emission from NGC 2623 is character- 
ized by a hard spectrum, suggesting that the scattered 
component of X-ray emission from a Compton thick (Njj 
> 10 24 cm" 2 ) AGN dominates the observed 2-10 keV 
X-ray flux (Maiolino et al. 2003). The intrinsic AGN lu- 
minosity is highly dependent on the unknown value of 
the scattering efficiency, and can be high if the efficiency 
is low. The optical non-detection of Seyfert signatures in 
NGC 2623 suggests that the putative AGN is a buried 
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one, whose scattering efficiency is expected to be low 
(Fabian et al. 2002). Evans et al. (2008) performed multi- 
wavelength observations of NGC 2623, and detected the 
high-excitation forbidden-emission line [NeV]14.3/xm in 
the infrared, indicative of narrow line region clouds pho- 
toionized by AGN radiation. The intrinsic AGN lumi- 
nosity can be high because the narrow line regions are 
expected to be under-developed in a buried AGN (Iman- 
ishi et al. 2006a, 2007a). The high CO (1-0) flux at the 
core of NGC 2623 (Bryant & Scoville 1999) suggests that 
detection of nuclear HCN(l-O) emission is feasible. 

Mrk 266 (NGC 5256; z = 0.028) is a LIRG (L IR = 
1O 115 L0), consisting of two main merging nuclei sepa- 
rated by — 10 arcsec, a south-western (SW) and north- 
eastern (NE) nucleus (Mazzarella & Boroson 1993). Op- 
tical spectroscopy by various groups suggested that the 
SW nucleus shows clear Seyfert 2 signatures but the NE 
one does not (Osterbrock & Dahari 1983; Kollatschny 
& Fricke 1984; Mazzarella & Boroson 1993; Wu et al. 
1998; Goncalves et al. 1999; Ishigaki et al. 2000). How- 
ever, the most systematic optical spectral classification 
of galaxies classified the SW and NE nuclei of Mrk 266 as 
LINER and Seyfert 2, respectively (Veillcux et al. 1995). 
High spatial resolution X-ray data obtained with Chan- 
dra found that the X-ray spectrum of the NE nucleus 
is characterized by a heavily-obscured AGN, while that 
of the SW nucleus consists of an obscured AGN and a 
strong starburst component (Brassington et al. 2007). 
The detection of a strong 3.3-/im PAH emission feature 
(Imanishi 2002), a good starburst indicator (Moorwood 
1986; Imanishi & Dudley 2000), supports the presence of 
a strong starburst in the SW nucleus. The radio 8.44- 
GHz (3.55 cm) to infrared (1-1000 /im) luminosity ratio 
of Mrk 266 is a factor of 2-3 higher than other LIRGs, 
suggesting the presence of a radio- intermediate (or radio- 
loud) AGN (Condon et al. 1991). The radio emission 
of the NE nucleus is dominated by a spatially compact 
(<0'.'3) component and is much brighter than that of the 
SW nucleus (Condon et al. 1991). The infrared 5-35 /im 
spectrum, taken with Spitzer IRS, is dominated by large 
equivalent-width PAH emission, which is typical of star- 
burst galaxies (Brandl et al. 2006). However, we note 
that the SW and NE nuclei arc not sufficiently covered 
with SL (5.2-14.5 /im) and LL (14-38 /j,m) spectra, re- 
spectively (Brandl et al. 2006), possibly missing AGN 
signatures, given the lack of full Spitzer spectral cover- 
age for both nuclei. 

Arp 193 (z=0.023) is a LIRG (L m = 1O 116 L ) with 
a long narrow tail along the south-eastern direction 
(Rothberg & Joseph 2004). It is classified optically 
as a LINER (Veilleux et al. 1995), showing no obvi- 
ous optical AGN signatures. The mid-infrared 8-20 
/im dust emission is dominated by a spatially-extended 
component with no prominent compact core (Soifer et 
al. 2001). The emission surface brightness is —2 x 
1O 12 L (Soifer et al. 2001), which is within the range 
of starburst phenomena. The HCN(l-O) and HCO+(l 
0) fluxes of Arp 193, measured with single-dish radio 
telescopes, are high (Solomon et al. 1992; Gracia-Carpio 
et al. 2008). Millimeter interferometric observations re- 
vealed that the CO(1-0) emission shows a strong nuclear 
compact component (Downes & Solomon 1998; Bryant 
& Scoville 1999). Hence, we expect that detection of 
nuclear HCN(l-O) emission is feasible if the nuclear to 



extended flux ratio for the HCN(l-O) line is similar to 
or higher than that of the CO (1-0) line. Papadopoulos 
(2007) found that the HCN(4-3) to HCN(l-O) flux ratio 
of Arp 193 is much lower than that of other LIRGs. The 
weak high-excitation HCN(4-3) line may suggest that 
the dominant energy source is a spatially-extended dif- 
fuse one, rather than a compact high-emission surface- 
brightness one. So far, no obvious AGN signatures have 
been seen in the Arp 193 observational data (Vega et al. 
2008). 

NGC 1377 (z=0.006) has an infrared luminosity of L m 
= 10 10 ^L©, and so is not a LIRG in a strict sense. This 
galaxy shows a deficit of 21-cm emission relative to in- 
frared emission (Roussel et al. 2003). Although it is op- 
tically unclassified (Veilleux et al. 1995), and thus no 
AGN signatures are present in the optical spectrum, the 
infrared 3 25 /im spectrum is dominated by a PAH-frec 
continuum (Imanishi 2006; Roussel et al. 2006), typi- 
cal of AGN-dominated galaxies. The 9.7-/im and 18- 
/im silicate dust absorption features are extremely strong 
(Roussel et al. 2006), suggesting that the putative AGN 
is deeply buried in dust. NGC 1377 is included in our 
sample to investigate whether even a buried AGN can- 
didate with low absolute infrared luminosity can have a 
high HCN(1-0)/HCO+(1-0) brightness-temperature ra- 
tio. This inclusion is important to test whether the 
general trend of high HCN(1-0)/HCO+(1-0) brightness- 
temperature ratios found in ULIRGs with buried AGN 
signatures (Imanishi et al. 2006b, 2007b) is indeed due 
to AGN-related phenomena, or simply the result of ex- 
treme molecular gas properties in ULIRGs (see §1). The 
CO(1-0) emission is detected with single-dish radio tele- 
scope observations (Roussel et al. 2003). If a significant 
fraction of the CO (1-0) emission comes from the nu- 
clear regions, and unless the nuclear HCN(1-0)/CO(1-0) 
brightness-temperature ratio is extremely low, then we 
estimate that detection of nuclear HCN(l-O) emission is 
feasible. 

3. OBSERVATIONS AND DATA REDUCTION 
3.1. Millimeter Inter -ferometry 

We performed interferometric observations of LIRGs, 
at HCN(l-O) (A rost = 3.3848 mm and v rcst = 88.632 
GHz) and HCO+(1-0) (A rost = 3.3637 mm or v rcst = 
89.188 GHz) lines, using the Nobeyama Millimeter Ar- 
ray (NMA) at the Nobeyama Radio Observatory (NRO) 
between 2006 December and 2008 March. For Mrk 266 
and NGC 1377, since no interferometric CO(1-0) maps 
had been published, we conducted NMA observations at 
CO(1-0) (A rost = 2.6026 mm or u rest = 115.271 GHz). 
Table 2 summarizes the detailed observing log. The 
NMA consists of six 10-m antennas and observations 
were undertaken using the AB (the longest baseline was 
351 m), C (163 m), and D (82 m) configurations. 

The backend was the Ultra- Wide-Band Correlator 
(UWBC) (Okumura et al. 2000) which can cover 1024 
MHz with 128 channels at 8-MHz resolution. For 
HCN(l-O) and HCO+(1-0) observations, the central fre- 
quency for each source (Table 2) was set to cover both the 
redshifted HCN(l-O) and HCO+(1-0) lines simultane- 
ously. A bandwidth of 1024 MHz corresponds to ~3500 
km s" 1 for the redshifted HCN(l-O) and HCO+(1-0) 
lines at v ~ 86-89 GHz and -2700 km s" 1 for the red- 
shifted CO(1-0) lines at v ~ 112-115 GHz. The fields 
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of view at 86-89 GHz and 112-115 HGz arc -77" and 
~62" at full- width at half- maximum (FWHM), respec- 
tively. The Hanning window function was applied to re- 
duce side lobes in the spectra, so that the effective reso- 
lution was widened to 16 MHz or 54 km s _1 (42 km s _1 ) 
at v ~ 86-89 GHz (112-115 GHz). 

The UVPROC-II package developed at NRO (Tsut- 
sumi et al. 1997) and the AIPS package of the National 
Radio Astronomy Observatory were used for standard 
data reduction. Corrections for the antenna baselines, 
band-pass properties, and time variation in the visibility 
amplitude and phase were applied to all of the data (Ta- 
ble 2). After discarding a fraction of the data with large 
phase scatter due to poor millimeter seeing and clipping 
a small fraction of data with unusually high amplitude, 
the data were Fourier-transformed using a natural UV 
weighting. The flux calibration was made using observa- 
tions of Uranus or appropriate quasars whose flux levels 
had been measured at least every month in the NMA 
observing seasons (Table 2). A conventional CLEAN 
method was applied to deconvolve the synthesized beam 
pattern. Table 3 summarizes the total net on-source in- 
tegration times and synthesized beam patterns. 

3.2. Ground-based Infrared K- (1.9-2.5 fim) and L- 
(2.8-4-1 fim) band Spectroscopy 

To put stronger constraints on the hidden energy 
sources at these millimeter-observed LIRG nuclei, we 
performed ground-based infrared K- (1.9-2.5 /im) and L- 
(2.8-4.1 /im) band spectroscopy of NGC 2623, Mrk 266, 
and Arp 193, as their L-band spectra were not available 
in the literature. NGC 1377 was not observed, because 
its L-band spectrum was given by Imanishi (2006). 

In short, through infrared L-band spectroscopy, we can 
investigate the energy sources based on the equivalent 
width of the 3.3 /an PAH emission and the optical depths 
of absorption features at A rcs t ~ 3.05 /jm by ice-covered 
dust g 3.4 /im by bare carbonaceous 

dust grains (Imanishi & Dudley 2000; Imanishi & Mal- 
oney 2003; Imanishi et al. 2006a). A normal starburst 
galaxy should always show large equivalent- width 3.3 
/im PAH emission, while a pure AGN produces a PAH- 
free continuum (Imanishi & Dudley 2000; Imanishi et 
al. 2006a). For absorption features, optical depths have 
upper limits in a normal starburst in which stellar en- 
ergy sources and dust are spatially well mixed, while the 
depths can be arbitrarily large in a buried AGN with 
a more centrally concentrated energy source geometry 
than the surrounding dust (Imanishi & Maloncy 2003; 
Imanishi et al. 2006a). In the infrared if-band spectra, 
emission from stars and AGNs is distinguishable using 
CO absorption features at A ros t = 2.3-2.4 /im because 
the CO absorption features are produced by stars older 
than 10 6 yr, but not by AGN-heated hot dust emission. 

Infrared K- (1.9-2.5 /mi) and L- (2.8-4.1 /mi) band 
spectra were taken using SpeX (Rayner et al. 2003) at- 
tached to the IRTF 3-m telescope atop Mauna Kea, 
Hawaii, on 2008 Apr 19 and 20 (UT). The 1.9-4.2 /im 
cross-dispersed mode was employed, so that K- and L- 
band spectra were obtained simultaneously. We chose 
a narrow O'.'8-wide slit to pinpoint the infrared contin- 
uum emission peaks of individual LIRG nuclei where the 
putative AGN may be located. The resulting spectral 
resolution using this slit is R — 1000 in the K- and L- 



bands. 

Photometric conditions persisted throughout the ob- 
servations. The seeing at K was measured in the range 
C/6-078 in FWHM. A standard telescope nodding tech- 
nique (ABBA pattern) with a throw of 7'/5 was employed 
along the slit. Each exposure was 15 sec, and two coadds 
were made at each position. The telescope tracking was 
monitored with the SpeX infrared slit-viewer. Table 4 
summarizes the detailed observing log. 

Appropriate standard stars (Table 4), with an airmass 
difference of <0.1 for individual LIRGs, were observed 
to correct for the wavelength-dependent transmission of 
Earth's atmosphere. The K- and i-band magnitudes of 
the standard stars were estimated based on their U-band 
(0.6 /xm) magnitudes, and V — K and V — L colors, of 
the corresponding stellar types (Tokunaga 2000). The 
estimated if -band magnitudes agree very well with the 
2MASS measurements (Skrutskie et al. 2006). 

Standard data reduction procedures were employed us- 
ing IRAF 3 . Initially, frames taken with an A (or B) beam 
were subtracted from frames subsequently taken with a 
B (or A) beam, and the resulting subtracted frames were 
added and divided by a spectroscopic flat image. Then, 
bad pixels and pixels hit by cosmic rays were replaced 
with interpolated values from surrounding pixels. Fi- 
nally the spectra of LIRG nuclei and standard stars were 
extracted by integrating signals over l'/8-2'/7, depend- 
ing on actual signal profiles. Wavelength calibration was 
performed based on standard star spectra which reflect 
the wavelength-dependent nature of Earth's atmospheric 
transmission curve with high signal-to-noise ratios. The 
LIRG nuclei spectra were divided by the observed spectra 
of standard stars, and then multiplied by the spectra of 
blackbodies with temperatures appropriate to individual 
standard stars (Table 4). 

Flux calibration was done based on signals of LIRGs 
and standard stars detected inside our slit spectra. Since 
we employed a narrow 0'.' 8- wide slit, our spectra are 
sensitive only to spatially-compact emission, and miss 
a spatially-extended component. To reduce data point 
scatter, appropriate binning of spectral elements was per- 
formed. 

4. RESULTS 

4.1. Millimeter Inter -ferometric Data 

For NGC 2623, Mrk 266, and Arp 193, our spectra at 
the HCN(l-O) or HCO+(1-0) emission peaks show that 
the flux levels between these lines are above zero, indi- 
cating that continuum emission is present. We combined 
data points that are unaffected by these lines and made 
interferometric maps of the continuum emission. Figure 
1 presents the contours of the continuum emission for 
NGC 2623, Mrk 266, and Arp 193, in the vicinity of the 
redshifted HCN(l-O) and HCO+(1-0) emission at 86-89 
GHz. For Arp 193, the continuum emission (~15 mJy) 
is clearly (>8er) detected. For NGC 2623, and Mrk 266 
SW and NE, signs of continuum emission are marginally 
seen (~3cr). Continuum emission was not detected in 
the HCN(1-0)/HCO+(1-0) data of NGC 1377, or the 

3 IRAF is distributed by the National Optical Astronomy Obser- 
vatories, operated by the Association of Universities for Research 
in Astronomy, Inc. (AURA), under cooperative agreement with 
the National Science Foundation. 
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CO(l-O) data of Mrk 266 and NGC 1377. 

Figure 2 displays integrated intensity maps of the 
HCN(l-O) and HCO+(1-0) emission of the observed four 
LIRGs. Figure 3 presents the interferometric maps of 
CO(l-O) emission for Mrk 266 and NGC 1377. Since 
LIRGs contain highly concentrated nuclear molecular gas 
(Sanders & Mirabel 1996), the fraction of high-density 
molecular gas (njj > 10 4 cm~ 3 ) is expected to increase 
(Gao & Solomon 2004). Such molecular gas is better 
probed with HCN(l-O) and HCO+(1-0) rather than the 
widely used CO(l-O) because the dipole moments of 
HCN and HCO + (/i > 3 dcbye) are much larger than 
CO (/i ~ 0.1 debyc; Botschwina et al. 1993; Millar ct 
al. 1997). Hence, the maps in Fig. 2 reflect the spatial 
distribution of the dense molecular gas phase. 

Mrk 266 has a double-nucleus morphology. The CO(l- 
0) emission is sufficiently bright to investigate its veloc- 
ity information for each individual nucleus, as shown in 
Figure 4. In this CO(l-O) channel map (Fig. 4), the 
NE nucleus shows a CO (1-0) emission peak at lower 
frequency (112.08-112.11 GHz) than the SW nucleus 
(112.17-112.19 GHz), suggesting that the NE nucleus is 
more redshifted than the SW nucleus. Previously ob- 
tained optical spectra also show that the NE nucleus is 
more redshifted than the SW nucleus (20-130 km s _1 ) 
(Mazzarella & Boroson 1993; Kollatschny & Fricke 1984; 
Kim ct al. 1995). 

Figure 5 shows spectra around the HCN(l-O) and 
HCO + (1-0) lines at the nuclear positions of the ob- 
served four LIRGs. Spectra around the CO (1-0) line 
for Mrk 266 SW and NE, and NGC 1377 are presented 
in Figure 6. For NGC 1377, HCN(l-O) emission might 
be marginally detected in the integrated intensity map 
close to the CO(1-0) emission peak (Figs. 2 and 3). In 
the spectrum at this peak position, five successive spec- 
tral elements at the expected HCN(l-O) frequency show 
flux excess, compared to the surrounding pixels (Fig. 5, 
lower- left). For reference, a spectrum with less binning, 
despite being slightly noisier, also displays signatures of 
HCN(l-O) emission (Fig. 5, lower-right). We thus regard 
the HCN(l-O) emission signature as real in NGC 1377. 

Figure 7 presents Gaussian fits to the detected HCN(1- 
0), HCO+(1-0), and CO(1-0) lines. In the HCN(1- 
0)/HCO + (1-0) spectra, the central velocity and line 
width of the Gaussian components are determined in- 
dependently for both lines. For NGC 1377, although 
HCN(l-O) and CO(1-0) emission are seen in similar po- 
sitions (Figs. 2 and 3), their velocity profiles are signifi- 
cantly different (Fig. 7). It may be that the HCN(l-O) 
line reflects high density molecular gas in the central nu- 
clear region of NGC 1377, while CO(1-0) emission probes 
the surrounding diffuse molecular gas toward the nuclear 
direction. We should take care interpreting the HCN(1- 
0)/CO(1-0) brightness-temperature ratio of NGC 1377 
because HCN(l-O) and CO (1-0) lines may probe physi- 
cally unrelated molecular gas components. 

Table 5 summarizes the derived Gaussian fitting pa- 
rameters. The derived fluxes of HCN(l-O) and HCO+(l- 
0) , and their brightness-temperature ratios at the nuclear 
peak positions (i.e., spatially unresolved component) are 
also summarized in Table 5. For NGC 2623 and Arp 
193, where continuum emission is subtracted, the inte- 
grated HCN(l-O) and HCO+(1-0) fluxes are also esti- 
mated from the peak contours of the integrated intensity 



maps in Fig. 2. The values estimated in this way agree 
with those based on the Gaussian fits to within 30%. Ta- 
ble 5 also includes CO(1-0) fluxes at the nuclear peaks, 
derived from our NMA data (Mrk 266 and NGC 1377) 
or taken from the literature (NGC 2623 and Arp 193). 

The CO (1-0) emission at Mrk 266 SW clearly displays 
spatially-extended structures (Fig. 3). The total CO(l- 
0) flux, including this spatially-extended component, is 
estimated to be 162 Jy km s _1 . 

CO(1-0) fluxes of Mrk 266 and NGC 1377, and 
HCN(l-O) and HCO+(1-0) fluxes of Arp 193, have been 
measured using single-dish radio telescopes (Sanders et 
al. 1986; Roussel et al. 2003; Gracia-Carpio et al. 2008). 
Table 6 compares these single-dish measurements with 
our NMA interferometric data. For Mrk 266 CO(1-0) 
emission, our interferometric flux measurements provide 
good agreement, within ~15%, to old single-dish mea- 
surements. For Arp 193, our HCN(l-O) flux is more com- 
parable to the latest IRAM 30m measurement (Gracia- 
Carpio et al. 2008), but is a factor of two smaller than 
the old IRAM one (Solomon et al. 1992), whose high 
flux measurement was noted by Gracia-Carpio et al. 
(2008). For Arpl93 HCO + (1-0) emission, our measure- 
ment agrees with the latest IRAM 30m data (Gracia- 
Carpio et al. 2008). For these lines, the bulk of the emis- 
sion should be covered in our interferometric data. For 
the nearest source, NGC 1377, however, our CO(1-0) 
flux measurement is more than a factor of four smaller 
than the SEST telescope measurement, suggesting that 
our interferometric data miss a spatially very extended 
component, and detect only nuclear compact molecular 
gas. 

4.2. Ground-based Infrared Spectra 

Figure 8 presents ground-based infrared K- (1.9-2.5 
/zm) and L- (2.8-4.1 /xm) band spectra of NGC 2623, 
Mrk 266 SW and NW, and Arp 193. if-band spectra 
with narrower wavelength coverage are available for NGC 
2623, Arp 193, and Mrk 266 SW (Ridgway et al. 1994; 
Goldader et al. 1997). For Mrk 266 SW, an i-band spec- 
trum is also available (Imanishi 2002). A if-band spec- 
trum of Mrk 266 NE, and L-band spectra of NGC 2623, 
Mrk 266 NE, and Arp 193 are first presented in this pa- 
per. For Mrk 266, this is the first infrared spectrum 
that clearly resolves emission from the Mrk 266 NE and 
SW nuclei, allowing separate investigation of the energy 
sources of individual nuclei for the first time. 

4.2.1. Continuum Flux Level 

For NGC 2623, and Mrk 266 NE and SW, the if-band 
flux levels of our SpeX spectra (078 slit width) are a 
factor of 1.7-2.7 smaller than the if-band photometric 
measurements with a 5" aperture (Carico et al. 1990). 
For NGC 2623, Mrk 266 SW, and Arp 193, the K-b&nd 
flux levels of our SpeX spectra (078 slit width) are 30- 
50% smaller than those measured in 3" x 9" aperture 
spectra (Goldader et al. 1997). The difference is great- 
est for Mrk 266 SW, whose spatially-extended if-band 
emission is strongest in our spectra along the slit direc- 
tion. The smaller flux levels in our spectra are reason- 
able, given that our smaller aperture misses a spatially 
extended emission component. 

4.2.2. Hydrogen Recombination Lines 
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Strong Paa (1.87 /im) and Bra (4.05 /zm) emission 
lines are visible in some spectra. Enlarged spectra 
around these strong emission lines are shown in Figure 9. 
In the if -band spectra, Br 7 (2.17 fim) and H 2 (l-0) S(l) 
(2.12 /zm) emission lines are detectable. Table 7 sum- 
marizes the strengths and line widths of these emission 
lines. If emission from the broad line regions close to the 
central SMBH in an AGN is strong and only modestly 
obscured, then we may be able to detect broad (FWHM 
> 1500 km s" 1 ) emission line components in the infrared, 
where dust extinction is much smaller than the optical, 
as was seen in optical Seyfert 2 galaxies (Hill et al. 1996; 
Veilleux et al. 1997; Lutz et al. 2002). However, we see 
no obvious signatures of broad line components. Since 
buried AGNs tend to contain a larger amount of obscur- 
ing dust than Seyfert 2 AGNs (Imanishi et al. 2006a, 
2007a, 2008), detection of the broad emission line com- 
ponents would be more difficult, even if the buried AGNs 
emit intrinsically luminous broad emission lines. 

In Table 7, since our small aperture spectra miss 
spatially-extended emission, the absolute flux is not 
meaningful. Only the flux ratios between different emis- 
sion lines measured with the same apertures will be used 
in our discussion. 

4.2.3. PAH Emission 

The 3.3-/im PAH emission is clearly seen in the L- 
band spectra of all the observed LIRGs. To estimate 
the 3.3-/xm PAH emission strength, we make the rea- 
sonable assumption that the profiles of the 3.3-/im PAH 
emission in these LIRGs are similar to those of Galactic 
star-forming regions and nearby starburst galaxies (type- 
1 sources in Tokunaga et al. 1991), following Imanishi ct 
al. (2006a). The adopted profile reproduces the observed 
3.3-/zm PAH emission features of the LIRGs reasonably 
well. Table 8 summarizes the fluxes, luminosities, and 
rest-frame equivalent widths of the 3.3 zxm PAH emission 
feature. Only the equivalent width is not significantly af- 
fected by flux loss in our small-aperture spectra, and so 
will be used in our discussion. 

4.2.4. CO Absorption 

Figure 10 shows enlarged spectra in the longer wave- 
length portion of the K-b&nd spectra. All spectra show 
spectral gaps in the continuum at A Q b s > 2.35 /im in 
the observed frame. We attribute the gaps to CO ab- 
sorption features at A rost = 2.31-2.4 /im produced by 
stars older than 10 6 yr (Oliva et al. 1999; Ivanov ct 
al. 2000; Imanishi & Alonso-Hcrrcro 2004; Imanishi & 
Wada 2004). To estimate the CO absorption strengths, 
we adopt the spectroscopic CO index (CO spoc ) defined 
by Doyon et al. (1994) and follow the procedures previ- 
ously applied to other LIRGs by ourselves (Imanishi ct 
al. 2004; Imanishi & Nakanishi 2006). Power-law contin- 
uum levels (Fa = a x X 13 ), shown as solid lines in Fig. 
10, are determined using data points at A rcst = 2.05- 
2.29 /im, excluding obvious emission lines. The derived 
COgpcc values are summarized in column 5 of Table 8. 

Ridgway et al. (1994) also estimated CO spoc = 
0.24±0.03 for NGC 2623, with a 2'.'7 aperture, which 
agrees with our measurement of CO spC c = 0.27 (Table 8) 
to within uncertainty. Goldader et al. (1997) measured 
the photometric CO index inside 3" x 9" apertures and 



estimated CO ph = 0.23, 0.16, and 0.18 for NGC 2623, 
Mrk 266SW, and Arp 193, respectively. Using the for- 
mula CO spoc = 1-46 x COp h - 0.02 (Goldader et al. 
1997), these CO p hot values are converted to CO spoc = 
0.32, 0.21, and 0.24 for NGC 2623, Mrk 266SW, and Arp 
193, respectively, which are slightly (0.02-0.05) larger 
than our measurements (Table 8). 

5. DISCUSSION 

5.1. Comparison of HCN(1-0)/HCO+ (1-0) 
Brightness-Temperature Ratios with Other Galaxies 

Figure 11 plots the HCN(1-0)/HCO+(1-0) and 
HCN(1-0)/CO(1-0) brightness-temperature ratios for 
the four LIRGs. Previously obtained data points of 
nearby LIRGs (Imanishi et al. 2004, 2006b; Imanishi & 
Nakanishi 2006; Imanishi et al. 2007b), starbursts, and 
Seyfert galaxies (Kohno 2005) are also plotted. As stated 
by Imanishi ct al. (2006b), the HCN(1-0)/HCO+(1-0) 
brightness-temperature ratios in the ordinate are mainly 
used in our discussions for the following reasons. First, 
since HCN(l-O) and HCO+(1-0) have similarly high 
dipole moments (it > 3 debye; Botschwina et al. 1993; 
Millar et al. 1997), it is very likely that similar high- 
density molecular gas is probed with both lines. Next, 
as both the HCN(l-O) and HCO+(1-0) lines are observed 
simultaneously with the same NMA array configuration, 
their beam patterns are virtually identical. Finally, both 
HCN(l-O) and HCO+(1-0) fluxes are measured' at the 
same time with the same receiver and same correlator 
unit, under the same weather conditions, so that possible 
absolute flux calibration uncertainties of intcrferometric 
data do not propagate to the ratio, which is dominated 
by statistical noise and fitting errors (see Fig. 7). 

In Fig. 11, we note the following points: first, NGC 
1377 exhibits a high HCN(1-0)/HCO+(1-0) brightness- 
temperature ratio, as seen in AGN-dominated galaxy nu- 
clei. Second, the Mrk 266 NE nucleus shows a signif- 
icantly higher ratio than the SW nucleus. Finally, the 
HCN(1-0)/HCO+(1-0) brightness-temperature ratios of 
Mrk 266 SW and Arp 193 are low, as observed in star- 
burst galaxies, and that of NGC 2623 is higher than Mrk 
266 SW and Arp 193. 

5.2. Buried AGN Signatures in Individual Objects 

In this subsection, we look for buried AGN signatures 
in the infrared spectra of individual LIRG nuclei and 
then investigate their HCN(1-0)/HCO+(1-0) brightness- 
temperature ratios. In doing so, we have to account 
for the fact that starburst activity surrounds the cen- 
tral compact buried AGN, so starburst emission is less 
obscured by gas and dust than the buried AGN emission. 
The contribution from the buried AGN to the observed 
infrared flux can be small, even if the buried AGN is in- 
trinsically luminous and is energetically important, due 
to dust extinction. Thus, we need to examine the in- 
frared spectra carefully. 

5.2.1. Mrk 266 NE 

Although optical spectroscopy by various groups has 
failed to find AGN signatures in the Mrk 266 NE nucleus 
(Osterbrock & Dahari 1983; Kollatschny & Fricke 1984; 
Mazzarella & Boroson 1993; Wu et al. 1998; Goncalves et 
al. 1999; Ishigaki et al. 2000), buried AGN signatures are 
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discernible in our infrared spectra. First, Mrk 266 NE 
shows a high H 2 (l-0) S(l) to Br-f flux ratio, > f (Ta- 
ble 7). Observationally, starburst galaxies show H 2 (l-0) 
S(l) to Br 7 flux ratios of <f (Mouri & Taniguchi 1992), 
while high H 2 (l-0) S(l) to Brj flux ratios (>1) are found 
in optical AGNs (— Seyfert galaxies) (Mouri & Taniguchi 
1992), as well as buried AGN candidates (Imanishi et al. 
2004). Next, the rest-frame equivalent width of the 3.3- 
/im PAH emission of Mrk 266 NE (EW 3 . 3 pah ~ 30 nm; 
Table 8) is significantly smaller than that observed in 
starburst galaxies (~100 nm; Moorwood 1986; Imanishi 
& Dudley 2000), suggesting that a PAH-free AGN- heated 
hot dust continuum contributes significantly to the ob- 
served 3-4 /zm flux of Mrk 266 NE. 

The HCN(1-0)/HCO+(1-0) brightness-temperature 
ratio of Mrk 266 NE is also higher than that of the Mrk 
266 SW nucleus. A putative buried AGN, suggested from 
X-ray data (Brassington et al. 2007), may enhance the 
HCN(l-O) emission in Mrk 266 NE. 

5.2.2. NGC 1311 

The HCN(1-0)/HCO+(1-0) brightness-temperature 
ratio of NGC 1377 is high, as found in AGNs. The in- 
frared L-band spectrum is also dominated by PAH-free 
continuum emission (Imanishi 2006), as usually seen in 
AGNs. The strong 9.7 /im silicate dust absorption fea- 
ture detected in the Spitzer IRS infrared 5-35 ^.m spec- 
trum (Roussel et al. 2006) is also incompatible with a 
normal starburst, where stellar energy sources and dust 
are spatially well mixed. It requires a buried AGN- type 
centrally concentrated energy source geometry (Imanishi 
et al. 2007a). These overall observational results are nat- 
urally explained by the presence of a luminous buried 
AGN in NGC 1377. 

A PAH-free continuum and strong 9.7-/^,m silicate dust 
absorption could be explained by an exceptionally cen- 
trally concentrated extreme starburst whose emitting 
volume is predominantly occupied with HH-rcgions, with 
virtually no molecular gas and photo-dissociation regions 
(Fig. le of Imanishi et al. 2007a). Unlike ULIRGs with 
Lir > 10 12 L Q , the absolute infrared luminosity of NGC 
1377 is only ^1O 1O1 L , which could be accounted for by 
a small number of super star cluster whose emitting size 
is very small (<<100 pc) and emission surface bright- 
ness is high (Gorjian et al. 2001). In fact, Roussel et 
al. (2003) preferred such a compact starburst (super star 
cluster) scenario, but ruled out the possibility of an ener- 
getically important buried AGN, based on the following 
arguments: (1) estimated SMBH mass is too small to 
account for the luminosity of NGC 1377 with AGN ac- 
tivity, and (2) high excitation forbidden emission lines, 
usually seen in Seyfert galaxies, are not detected. 

Regarding the first point, the SMBH mass was esti- 
mated to be <2 x 10 5 M Q , from the observed radio 20- 
cm (1.5-GHz) continuum flux and the small line width 
of the infrared H 2 (l-0) S(l) emission line (FWHM < 25 
km s _1 ) (Roussel et al. 2003). The measured line width 
of millimeter CO (1-0) lines in our NMA spectrum (Fig- 
ure 6) is ^100 km s _1 in FWHM. After correction for 
the velocity resolution of NMA data (~42 km s _1 ; §3.1), 
we obtain the intrinsic CO(1-0) line width of ^90 km 
s- 1 in FWHM, which is similar to the SEST 15-m tele- 
scope measurement (Roussel et al. 2003). Assuming the 
CO (1-0) line width and SMBH mass relation given by 



Shields et al. (2006), the measured CO(1-0) line width 
provides a similar SMBH mass of -2 x 10 5 M Q . The Ed- 
dington luminosity of this SMBH mass is ^7 x 10 9 L Q . 
This is only slightly below the observed infrared lumi- 
nosity (~1 x 10 10 L Q ). In a buried AGN, the infrared 
luminosity should be comparable to the bolometric lu- 
minosity, because almost all of the energetic radiation 
from the central AGN is absorbed by the surrounding 
dust and re-emitted as infrared thermal dust emission. 
Thus, even if current constraints on the supermassive 
black hole mass are accurate, it is still possible that the 
bulk of the observed luminosity in NGC 1377 comes from 
buried AGN activity. 

However, uncertainties in the estimated SMBH mass 
in NGC 1377 could exist. First, the observed 20-cm flux 
from a buried AGN can be severely attenuated by free- 
free absorption, possibly leading to an underestimate of 
the SMBH mass. Second, the line width of H 2 (l-0) S(l) 
emission is not well calibrated for an SMBH mass esti- 
mate. Finally, in the small line width range of millimeter 
CO (1-0) (<100 km s" 1 ), the SMBH mass estimated from 
CO (1-0) line width yields systematically smaller results 
than other measurements (Shields et al. 2006). Hence, 
the actual SMBH mass in NGC 1377 could be higher 
than above estimates. 

Regarding the second point, the non-detection of high 
excitation emission lines is a natural consequence of a 
buried AGN because the central AGN is obscured by 
dust and gas along virtually all directions at the inner 
part (<1 pc), producing virtually no narrow line regions 
(= the main emitting sources of high excitation lines). 

Therefore, none of the currently available observational 
results preclude the presence of a luminous buried AGN 
in NGC 1377. Strong H 2 emission is observed in NGC 
1377 (Roussel et al. 2006). In general, H 2 emission 
is stronger in AGNs than starburst galaxies (Mouri & 
Taniguchi 1992). Although a major galaxy merger could 
produce strong H 2 emission (Van der Werf et al. 1993), 
NGC 1377 shows no sign of a major merger (Roussel et 
al. 2006). Energy source obscuration for NGC 1377 is ex- 
tremely high (Roussel et al. 2006), and the observed low 
radio to infrared luminosity ratio (Roussel et al. 2003) 
could be explained by severe flux attenuation of the ra- 
dio 20-cm (1.5 GHz) emission by free-free absorption. 

5.2.3. NGC 2623 

The EW3.3PAH and CO spcc values, and the H 2 (l-0) 
S(l) to Br7 flux ratio of NGC 2623 are all within star- 
burst range. The HCN(1-0)/HCO+(1-0) brightness- 
temperature ratios are slightly higher than Arp 193 and 
Mrk 266 SW, which may be due to an HCN(l-O) emis- 
sion enhancement by an X-ray detected buried AGN (see 
§2). 

5.2.4. Mrk 266 SW and Arp 193 

The large EW3.3PAH and CO spcc values, the small 
H 2 (l-0) S(l) to Br7 flux ratios, and the small HCN(1- 
0) /HCO + (1-0) brightness-temperature ratios of Mrk 266 
SW and Arp 193 are all explained by starburst activity 
only, with no significant AGN contribution required. 

5.3. Interpretation of High HCN(1-0)/HCO+ (1-0) 
Brightness-Temperature Ratios in Buried AGN 
Candidates 
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One natural explanation for the high HCN(1- 
0)/HCO + (1-0) brightness-temperature ratios found in 
AGNs is an HCN abundance enhancement. If molec- 
ular gas consists of small dense gas clumps with low 
volume filling factor, as widely supported from obser- 
vations (Solomon et al. 1987), an enhanced HCN abun- 
dance can result in higher HCN(l-O) flux regardless of 
whether the HCN (1-0) emission is optically thin or thick 
(Imanishi et al. 2007b). Several chemical calculations of 
HCN abundance in molecular gas around UV- and X-ray 
emitting energy sources have been published (Meijerink 
et al. 2006; Lintott & Viti 2006). Although an HCN 
abundance enhancement around an X-ray emitting en- 
ergy source (i.e., AGN) is predicted in some parameter 
ranges that are realistic for molecular gas around buried 
AGNs in LIRGs (high-F x range of Table 3 in Meijerink 
et al. 2006), an HCN abundance decrease is suggested 
in other reasonable parameter ranges (Meijerink et al. 
2006). The chemical calculation results are highly depen- 
dent on parameters, and it is currently unclear whether 
an HCN abundance is indeed enhanced in molecular gas 
in the close vicinity of buried AGNs in LIRG nuclei. 

The second possible explanation for the strong HCN(1- 
0) emission in AGNs is an infrared radiative pumping sce- 
nario (Aalto et al. 1995; Garcia-Burillo et al. 2006; Weiss 
et al. 2007). It is usually assumed that molecular gas is 
excited by collision. However, the HCN molecule has a 
line at infrared 14 /im that can be excited by absorbing 
photons at A ~ 14 /im. HCN(l-O) emission in the mil- 
limeter wavelength range could be enhanced through a 
cascade process. Since the emission surface brightness of 
an AGN is high, much of the surrounding dust is heated 
to several 100 K, producing strong mid-infrared 10-20 
jjm emission (§1). Hence, this infrared radiative pump- 
ing scenario for HCN could work effectively in AGNs. 

HCO + also has a line at infrared 12 fim, so that the 
infrared radiative pumping mechanism could work in a 
similar way. Observationally, the 14 /im HCN absorption 
features are detected in highly obscured LIRGs (Lahuis 
et al. 2007), but the 12 fim. HCO + absorption features 
are not (Farrah et al. 2007). It may be that this infrared 
radiative pumping scenario is working more effectively 
for HCN than HCO+, possibly enhancing the HCN(1- 
0)/HCO + (1-0) brightness-temperature ratios in AGNs. 

To test this infrared radiative pumping scenario, in 
Figure 12 we compare the observed HCN(1-0)/HCO+(1- 
0) brightness-temperature ratio with the infrared emis- 
sion surface brightness estimated by Soifer et al. (2000, 
2001) and Evans et al. (2003). We use the infrared 
emission surface brightness, rather than infrared lumi- 
nosity (Gracia-Carpio et al. 2006), because the bulk of 
the observed normal LIRGs with <10 12 L© (i.e., non- 
ULIRGs) are chosen because they display luminous 
buried AGN signatures. Although Gracia-Carpio et al. 
(2008) found an increasing trend of HCN(1-0)/HCO+(1 
0) brightness-temperature ratios with increasing galaxy 
infrared luminosities, our heterogeneous sample selection 
of normal LIRGs (non-ULIRGs) could be biased to AGN 
candidates (= strong HCN emitters) and artificially elim- 
inate this trend. The infrared radiative pumping sce- 
nario can, in principle, enhance the HCN(1-0)/HCO + (1 - 



0) brightness-temperature ratio, even for a low absolute 
infrared luminosity galaxy (e.g. NGC 1377), if infrared 
emission surface brightness is high. 

In Fig. 12, we may see a weak correlation but the scat- 
ter is large and the total sample size is small. Based on 
our current dataset, we cannot determine whether the 
infrared radiative pumping scenario is indeed at work in 
buried AGNs at LIRG's nuclei. Further detailed theoret- 
ical calculations that realistically incorporate the actual 
energy level populations of HCN and HCO + , as well as 
the clumpy structure of molecular gas around an AGN 
(Yamada et al. 2007), are needed to properly interpret 
the origin of strong HCN(l-O) emission in AGNs. Com- 
bination with higher transition lines in the submillimeter 
wavelength range (Wilson et al. 2008) will also help bet- 
ter understand the physical properties of molecular gas 
in LIRG's nuclei. 

6. SUMMARY 

We presented the results of millimeter interferomctric 
simultaneous HCN (1-0) and HCO + (1-0) observations of 
four LIRGs using NMA. When combined with previously 
observed LIRGs, ours is the largest LIRG interferometric 
HCN(l-O) and HCO+(1-0) survey. From our interfero- 
metric data, we extracted the HCN(l-O) and HCO+(1-0) 
fluxes at the nuclei, where putative AGNs are expected 
to be present. We derived the HCN(1-0)/HCO+(1-0) 
brightness-temperature ratios of the observed LIRG nu- 
clei and compared them with the ratios found in AGNs 
and starburst galaxies. The main result of this paper is 
the confirmation of the previously discovered trend that 
LIRGs with luminous buried AGN signatures at other 
wavelengths tend to show higher HCN(1-0)/HCO+(1- 
0) brightness-temperature ratios than those without. 
This reinforces the utility of the HCN(1-0)/HCO+(1-0) 
brightness-temperature ratio as a powerful observational 
tool for discovering elusive AGNs deeply buried in dust 
and molecular gas. 
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TABLE 1 

Detailed information on the observed LIRGs 



Object 


Redshift 


fl2 


/as 


/eo 


/ioo 


log L m (log L m /L e ) 


Far-infrared 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


(ergs s -1 ) 


Color 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 2623 (Arp 243) 


0.018 


0.21 


1.74 


23.13 


27.88 


45.1 (11.5) 


0.08 (cool) 


Mrk 266 (NGC 5256) 


0.028 


0.23 


0.98 


7.34 


11.07 


45.0 (11.5) 


0.13 (cool) 


Arp 193 (IC 883, UGC 8387) 


0.023 


0.26 


1.36 


15.44 


25.18 


45.2 (11.6) 


0.09 (cool) 


NGC 1377 


0.006 


0.44 


1.81 


7.25 


5.74 


43.7 (10.1) 


0.25 (warm) 



Note. — Col.(l): Object name. Col. (2): Redshift. Cols.(3)-(6): f 12 , f25, feo, and fi 00 arc IRAS FSC fluxes at 12, 25, 60, and 100 ^m, 
respectively. Col. (7): Decimal logarithm of the infrared (8—1000 M m ) luminosity in ergs s _1 calculated as follows: Lir — 2.1 X 10 39 X D(Mpc) 2 X 
(13.48 x fi2 4- 5.16 x /2s + 2.58 x /eo + /100) ergs s~ (Sanders & Mirabel 1996). The values in parentheses arc the decimal logarithms of the 
infrared luminosities in units of solar luminosities. Col. (8): IRAS 25 fim to 60 /xm flux ratio (f25/feo)- LIRGs with f25/f60 < (>) 0.2 arc classified 
as cool (warm) (Sanders ct al. 1988b). 



TABLE 2 





Ob 


SERVING LOG OF NMA MILLIMETEF 


: INTERFEROMETR 


IC OBSERVATIONS 






Object 
(1) 


Array 
configuration 
(2) 


Observing Date 
(UT) 

(3) 


Central 
frequency 
(4) 


Phase 
calibrator 
(5) 


Band-pass 
calibrator 
(6) 


Flux 
calibrator 
(7) 


NGC 2623 


C 


2007 Dec 13, 14, 15, 16 


87.30 


0851+202 


3C 84 


3C 345 


Mrk 266 


c 

D 

C (CO) 


2008 Jan 8, 9, 10 

2008 Feb 20, 21, 22, Mar 4, 5 

2007 Dec 19, 23 


86.52 
112.43 


1418+546 
1418+546 


30 84 
3C 84 


1749+096 
1749+096 


Arp 193 


c 

D 


2007 Dec 20 

2008 Feb 18, 19 


86.89 


1308+326 


3C 84 


1749+096 


NGC 1377 


D 
AB 

C 

C (CO) 


2006 Dec 12, 13 

2007 Jan 27, 28 

2007 Dec 13, 14, 15, 16, 19, 20 
2007 Dec 10, 11 


88.47 
114.95 


0420-014 
0420-014 


3C 454.3 
3C 454.3 


Uranus 
Uranus 



NOTE. — Col.(l): Object name. Col. (2): NMA array configuration. The mark "(CO)" means CO(1-0) observations, which were executed only 
for Mrk 266 and NGC 1377. Col. (3): Observing date in UT. Col. (4): Central frequency used for the observations. Col. (5): Object name used as a 
phase calibrator. Col. (6): Object name used as a band-pass calibrator. Col. (7): Object name used as a flux calibrator. 



TABLE 3 
Parameters of final NMA maps 



Object 

(i) 



Line 
(2) 



On source 
integration (hr) 
(3) 



Beam size 
(arcscc X arcscc) 
(4) 



Position angle of 
the beam (°) 
(5) 



NGC 2623 
Mrk 266 

Arp 193 
NGC 1377 



HCN/HCO+ 
HCN/HCO+ 
CO 

HCN/HCO+ 
HCN/HCO+ 
CO 



10 

34 
5 
8 

23 
4 



5.4 X 4.2 

4.3 x 3.4 

4.3 x 3.4 

7.1 x 6.6 

7.1 x 4.0 

6.7 x 3.6 



-48.4 
-31.5 
-31.5 
-52.5 
-13.6 
-18.8 



Note. — Col.(l): Object name. Col. (2): Observed line. HCN(1-0)/HCO+ (1-0) or CO(1-0). Col. (3): Net on-source integration time in 
hours. Col. (4): Beam size in arcsec x arcsec. Col. (5): Position angle of the beam pattern. It is 0° for the north-south direction, and increases 
counterclockwise on the sky plane. 
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TABLE 4 

Observing log of infrared spectroscopy with IRTF SpeX 



Object 


Date 


Integration 




Standard Stars 






(UT) 


(Min) 


Name 


K-mag 


L-mag 


Type 


T cS (K) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 2623 


2008 April 20 


60 


HR 3262 


3.9 


3.9 


F6V 


6400 


Mrk 266 SW 


2008 April 19 


40 


HR 4767 


4.9 


4.8 


F9V-G0V 


6000 


Mrk 266 NE 


2008 April 19 


60 


HR 4767 


4.9 


4.8 


F9V-G0V 


6000 


Arp 193 


2008 April 20 


60 


HR 4845 


4.5 


4.5 


GOV 


5930 



Note. — Column (1): Object name. Col. (2): Observing date in UT. Col. (3): Net on-source integration time in min. Col. (4): 
Standard star name. Col. (5): Adopted K-band magnitude. Col. (6): Adopted L-band magnitude. Col. (7): Stellar spectral type. 
Col. (8): Effective temperature. 



TABLE 5 

Gaussian fitting parameters and nuclear fluxes of HCN(l-O), HCO+(1-0), and CO(1-0) emission lines 



Object 




LSR velocity 






Flux 




HCN(1-0)/HCO+(1-0) 






[km s _1 ] 






[Jy km s _1 ] 




Ratio 




HCN(l-O) 


HCO+(1-0) 


CO(1-0) 


HCN(l-O) 


HCO+(1-0) 


CO(1-0) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 2623 


5570 


5620 




5.8 


4.0 


153 


1.5 


Mrk 266 SW 


8290 


8290 


8370 


1.0 


3.5 


88 


0.3 


Mrk 266 NE 


8330 


8370 


8250 + 8501 


1.1 


0.8 


18 


1.4 


Arp 193 


6905 + 7110 


6850 + 7130 




10.8 


12.7 


202 


0.9 


NGC 1377 


1500 




1740 


2.2 


<1.2 


9.6 


>1.8 



Note. — Col.(l): Object name. Cols.(2)-(4): LSR velocity {v opt = - 1) x c} of the HCN(l-O), HCO+(1-0), and CO(1-0) in [km s -1 ]. 
Cols.(5)-(7): Flux of the HCN(l-O), HCO+(1-0), and CO(1-0) in [Jy km s -1 ] at the nuclear peak position. For NGC 2623 and Arp 193, the 
CO(1-0) fluxes are adopted from Bryant & Scoville (1999). Col. (8): HCN(1-0)/HCO+ (1-0) ratio in brightness temperature (ex A 2 X flux density). 
The ratio is not affected by possible absolute flux calibration uncertainties in the NMA data (see §5.1). 



TABLE 6 

Comparison of HCN(l-O), HCO+(1-0), and CO(1-0) flux between our NMA interferometric measurements and single-dish 

telescope's measurements in the literature 



Nucleus 


Line 


Flux 


Reference 


(1) 


(2) 


(3) 


(4) 


Mrk 266 


CO(l-O) 


180 (SW + NE) a 


This work 






206 b 


Sanders ct al. (1986) 


Arp 193 


HCN(l-O) 


10.8 


This work 




25 c 


Solomon ct al. (1992) 






6.3 d 


Gracia-Carpio et al. (2008) 




HCO+(1-0) 


12.7 


This work 






10.0 d 


Gracia-Carpio ct al. (2008) 


NGC 1377 


CO(l-O) 


9.6 


This work 




47 ° 


Roussel et al. (2003) 



Note. — Col.(l): Object name. Col.(2): HCN(l-O), HCO+(1-0), or CO(1-0) line. Col.(3): Flux in [Jy km s -1 ]. Col.(4): 
Reference. 

a For the CO(1-0) flux from the SW nucleus, we adopt the total flux of 162 Jy km s _1 (§4.1), as the spatially-extended 
component is very strong. 

b A conversion factor of 42 Jy K _1 is assumed for FCRAO 14-m telescope (Kenny & Young 1988). 
C A conversion factor of 4.4 Jy K~ is assumed for IRAM 30-m telescope (Solomon et al. 1992). 

d A conversion factor of 6.0 Jy K _1 ("K" is antenna temperature in this case) is assumed for IRAM 30-m telescope (Gracia- 
Carpio et al. 2008). 

°A conversion factor of 27 Jy K _1 is assumed for SEST 15-m telescope at 115 GHz 
(http:/ /www.ls. eso.org/lasilla/Telescopes/SEST/html/telescope-instruments/telescope/index. html). 
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TABLE 7 

Properties of hydrogen emission lines 



Object 




Paa 




Bra 


Br7 


H 2 


Paa/Bra 


H 2 /Br 7 




flux 


FWHM 


flux 


FWHM 


flux 


flux 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 2623 


92 


515 


40 


510 


8 


6.5 


2.5 


0.8 


Mrk 266 SW 


40 


365 






7 


2.5 




0.3 


Mrk 266 NE 


50 


510 






2 


7 




3.5 


Arp 193 


72 


440 


39 


465 


11 


7.5 


2 


0.7 



Note. — Col.(l): Object name. Col.(2): Paa flux in [lCT 15 ergs s" 1 cm" 2 ]. Col. (3): Pact line width in FWHM in [km s" 1 ]. 
Col. (4): Bra flux in [1(T 15 ergs s" 1 cm" 2 ]. Col.(5): Bra line width in FWHM in [km s -1 ]. Col. (6): Br 7 flux in [TO -15 ergs s" 1 
cm" 2 ]. Col. (7): H 2 (l-0) S(l) flux in [HT 18 ergs s" 1 cm" 2 ]. Col. (8): Paa/Bra flux ratio. Col. (9): H 2 /Br7 flux ratio. 



TABLE 8 

Properties of 3.3-/*m PAH emission and 2.3-/im CO absorption 



Object 


I3.3PAH 


L3.3PAH 


EW3.3PAH 


COspcc 


(1) 


(2) 


(3) 


(4) 


(5) 


NGC 2623 


16 


1.0 


125 


0.27 


Mrk 266 SW 


9 


1.4 


70 


0.19 


Mrk 266 NE 


3 


0.4 


30 


0.24 


Arp 193 


21 


2.2 


140 


0.21 





Note. — Col.(l): Object name. Col. (2): Observed flux of 3.3-^m PAH emission in [10 14 ergs s 1 cm" 2 ]. Col. (3): Observed 
luminosity of 3.3-/im PAH emission in [10 41 ergs s -1 ]. Col. (4): Rest-frame equivalent width of the 3.3-/im PAH emission in 
[nm]. Col. (5): Spectroscopic CO index defined by Doyon et al. (1994). 
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NGC 2623 
continuum 




RIGHT ASCENSION (J2DDD) 



Arp 193 
continuum 





36.D 35.5 35.D 34.5 

RIGHT ASCENSION (J2DDD) 



Mrk 266 
continuum 



13 38 1S.D 185 



18.0 17.5 17.D 16.5 
RIGHT ASCENSION (J20DD) 



Fig. 1. — Continuum maps of NGC 2623, Mrk 266, and Arp 193 at v ~ 86-89 GHz. The crosses show the coordinates of main nuclei. The 
coordinates in J2000 arc (08 h 38 m 24.09 s , +25°45'16'.'5) for NGC 2623, (13 ,i 38 m 17.36 s , +48°16'32'.'5) for Mrk 266 SW, (13 h 38 m 17.84 s , 
+48°16'41'.'l) for Mrk 266 NE, and (13 h 20 m 35.39 s , +34°08'2T.'l) for Arp 193. The nuclear coordinates of NGC 2623 and Arp 193 are 
adopted from interferometric CO(l-O) maps by Bryant & Scoville (1999), after conversion from B1950 to J2000, using NED. That of Mrk 
266 is estimated from Chandra X-ray data by Brassington et al. (2007). The contours are 1.3 X 3 mjy beam -1 for NGC 2623, 0.8 X 3 
mjy beam -1 for Mrk 266, and 1.5 X (4, 5, 6, 7, 8) mjy beam -1 for Arp 193. The continuum fluxes at the peak positions arc ~4 mjy, 
~2.8 mjy, ~2.3 mjy, and ~15 mjy for NGC 2623, Mrk 266NE, Mrk 266SW, and Arp 193, respectively. Beam patterns are shown in the 
small squares at the lower-left corners. 
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NGC 1377 
HCO+ 




392 39.0 38.8 3B.6 
RIGHT ASCENSION (J2000) 



392 39.0 38 3 3B.6 
RIGHT ASCENSION (J2000) 



Fig. 2. — Left: HCN(1— 0) emission map. Right: HCO + (l— 0) emission map. The crosses show the coordinates of main nuclei. For NGC 
2623, Mrk 266, and Arp 193, the same coordinates as used in Figure 1 are shown. For NGC 1377, the coordinate in J2000 is (03 h 36 m 39.12 s , 
— 20°54'04"4), estimated from optical B-band data by Roussel et al. (2006). Continuum emission is subtracted for NGC 2623 and Arp 193, 
but not for Mrk 266 and NGC 1377 because continuum levels arc different for the NE and SW nuclei for Mrk 266 and continuum emission 
is not detected for NGC 1377 (<2 mjy). For NGC 2623, the contours are 0.66 X (4, 5, 6, 7) Jy km s" 1 for HCN, and 0.73 X (4, 5, 6, 7) 
Jy km s- 1 for HCO+. For Mrk 266, the contours arc 0.62 X (4, 5) Jy km s" 1 for HCN, and 0.55 X (3, 4, 5, 6, 7, 8) Jy km s" 1 for HCO+. 
For Arp 193, the contours are 1.0 X (6, 8, 10) Jy km s" 1 for HCN, and 1.0 X (6, 8, 10, 12) Jy km s" 1 for HCO+. For NGC 1317, the 
contours are 0.5 X 4 Jy km s _1 for HCN, and 0.5 X 3 Jy km s -1 for HCO+. 




Fig. 3. — CO(l— 0) emission for Mrk 266 and NGC 1377. No continuum subtraction was attempted, because the continuum emission is 
undetected (<5 mjy and <6 mjy for Mrk 266 and NGC 1377, respectively) and negligible compared to the strong CO(1-0) emission. The 
contours are 2.85 x (3, 5, 7, 9, 11, 13, 15, 17, 19, 21, 23, 25, 27, 29) Jy km s" 1 for Mrk 266, and 1.21 x (5, 6, 7, 8) Jy km s" 1 for NGC 
1377. 
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Fig. 4.— Channel map of CO(l-O) emission for Mrk 266. The contours arc 12 X (-3, 3, 5, 7, 9, 11, 13, 15) mjy beam -1 . The rms noise 
level is ~12 mjy beam -1 . 
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Fig. 9. — Enlarged spectra around Paa and Bra emission lines. The abscissa is the observed wavelength in /im, and the ordinate is flux 
F\ in 10~ 15 W m~ 2 fim^ 1 . Gaussian fits are overplotted as dotted lines. 
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Fig. 12. — HCN(1-0)/HCO" 1 " (1— 0) brightness-temperature ratio (ordinate) and infrared emission surface brightness (abscissa) for LIRGs 
with available infrared emission surface brightness data (Soifer et al. 2000, 2001; Evans et al. 2003). Only LIRGs whose HCN(1— 0) and 
HCO + (1-0) lines were simultaneously observed with NMA, and the prototypical starburst galaxy M82 (Kohno et al. 2001), are plotted. 
The LIRG sample includes IRAS 08572+3915, UGC 5101, Mrk 231, Mrk 273, Arp 220, IRAS 17208-0014, VV 114 E, Arp 193, NGC 
2623, Arp 299 A and C, and NGC 4418 (Imanishi et al. 2004, 2006b; Imanishi & Nakanishi 2006; Imanishi et al. 2007b). For Arp 220, the 
HCN(1-0)/HCO+ (1-0) brightness-temperature ratio is the average value of the eastern and western nuclei (Imanishi et al. 2007b). 



